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Measurements made on dental amalgam in tension indicate that amalgam exhibits three types of 
viscoelastic phenomena: (1) instantaneous elastic strain, (2) retarded elastic strain (transient creep), 
and (3) viscous strain (steady state creep). The combination of elastic plus retarded strain can be repre- 
sented by an equation of the form e = Acr J rB 2 & 1 where A and B are functions of time but not of the 
stress, a. The viscous strain rate can be represented by an equation of the form e v —Kcr m where 
K and m are constants of the material. By applying a nonlinear generalization of the Boltzmann super- 
position principle to a general equation describing the creep behavior of amalgam, the results of 
creep tests can be directly related to the results of stress-strain tests. 
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1 . Introduction 

There is no information in the literature on the 
viscoelastic behavior of dental amalgam except the 
limited amount that can be obtained from stress-strain 
relationships [1, 2]. 1 The stress-strain data indicate 
that amalgam exhibits a nonlinear relation between 
stress and strain over the entire range of stress in- 
vestigated: this has been found to be true in both 
tension and compression. It has been further noted 
that the shapes of the stress-strain curves vary with 
the strain rate indicating that the strain developed is 
not only dependent on the applied stress but also upon 
time of application of the stress. The results of tensile 
stress-strain investigation by Rodriguez and Dickson 
[1] indicated that dental amalgam might be a visco- 
elastic material. It was, therefore, the objective of 
this study (1) to make an exploratory investigation of 
the types of viscoelastic phenomena exhibited by 
dental amalgam in tension, (2) to describe the visco- 
elastic phenomena exhibited by amalgam in terms of 
available viscoelastic theory and (3) to make avail- 
able a practical example and method of application 
of viscoelastic theory to a dental material. 

2. Theory 

In general a strain-hardened material will ex- 
hibit one or more of three aspects of viscoelastic phe- 
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nomena: (1) elastic behavior, (2) viscous behavior 
(steady-state creep), and (3) retarded elastic behavior 
(transient creep). 

In a material which exhibits all three phenomena 
linearly related to the applied stress and in which the 
responses are additive (obey a superposition relation), 
the strain at any time may be represented by the equa- 
tion 



./ CT+cr \ X J'(T)(l-e- t l*)dT+— , 
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e is the creep strain, 

J () a is the instantaneous elastic response, 

\J'{t) (1 — e~ t,T )dr is the retarded elastic re- 
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is the viscous response, 
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a is the applied stress, 

J is the elastic compliance, the reciprocal 
elastic modulus, 

./'(r) is a continuous distribution of retarded elas- 
tic compliances as a function of the variable 
retardation time r of a continuous distribu- 
tion of responses 

t is the time after application of the stress, and 

7] is the coefficient of viscosity. 

The relationship between creep strain and stress 
may not be linear as assumed above. For example, a 
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material may obey the following equation: 
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where K{<j — a {) ) n 't = for all values of cr ^ cr and 
(3(y) is a continuous distribution of retarded 

elastic compliances as a function of the 
variable retardation time y of a continuous 
distribution of responses in the nonlinear 
second degree stress response of the 
material, 
<j is the applied stress level at which the 

material begins to exhibit a viscous flow, 
and 
K and m are constants at any given temperature. 
Where the creep strain is linearly related to the 
applied stress, the strain can be written as a product 
of a general creep compliance, which is a function of 
time only, and the applied stress as given by the fol- 
lowing equation: 



e=J(t)a, 



(3) 



where 



e is the creep strain developed, 

cr is the applied stress, and 

J (t) is the general creep compliance. 

For example, in the linear creep behavior as seen in 
equation (1), the general creep compliance J (t) is 
of the form: 
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In the case that creep strain is linearly related to 
the applied stress, the stress-strain behavior for a 
material may be related to the creep behavior by 
means of the linear superposition principle devel- 
oped by Boltzmann [3] as shown below: 



(5) 



or 



E(a)= J(T-0)da(0). 



E(T)= f ' J(T-O) ( -^P-d6, 
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E(T) 



dcr(O) 
dO 



is the applied stress related to the ex- 
perimental time in a prescribed ex- 
perimental functional relation, 
is the strain of a stress-strain relation 
measured at a time T subsequent to the 
variable time 0, and 

dO is the increment of stress from time 
to time (0 + dO). 



For example, applying the Boltzmann superposi- 
tion principle to the example of creep given in eq (1) 
gives the following relation between creep behavior 
and the strain measured in a stress-strain test [4]: 
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A discussion of the application of the Boltzmann 
superposition principle as applied to linear viscoelastic 
phenomena is given by Leaderman [5] . However, if 
a material exhibits a nonlinear relation in its retarded 
strain (transient creep) and the applied stress, as 
discussed earlier, the linear Boltzmann superposition 
principle cannot be applied to relate strain of a stress- 
strain test on the material to the creep strain. A gen- 
eralization of Boltzmann's superposition principle 
developed by Nakada [6] may be used to relate non- 
linear retarded elastic creep behavior to the stress- 
strain behavior of the material as given: 
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Now for example, consider the application of the 
nonlinear superposition principle of Nakada as applied 
to the case of retarded strain (transient creep) related 
to the applied stress by a second degree equation as 
shown earlier: 
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where ty(t) and O(0 could be represented by the 
equation shown below in accordance with the Voigt 
model r 

^(t)=\J f (T)(l-€-^)dr, 
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If one applies an approximation of Nakada's general 
formulation described by Leaderman, McCrackin, 
and Nakada [7], the equation representing the re- 
tarded elastic response takes the following form: 



E (<r(T)) = 
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Considering an example of a material which exhibits 
all three phenomena but in which the creep strain is 
a nonlinear function of the applied stress according 
to eq (2), then applying the above approximation of 
the nonlinear superposition principle, the stress- 
strain behavior of the material may be related to the 
creep behavior by the following equation: 
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However, as will be discussed later, it was found in 
the investigation that a different approximation is 
required to describe the behavior of dental amalgam. 



was taken immediately. Strain readings were then 
made at 15 s intervals for 4 min and at increasingly 
longer intervals until the strain rate became con- 
stant (usually after approximately 1.5 hr). At the end 
of this period, the load was removed, a strain reading 
was taken immediately and the recovery curve was 
followed by reading first at 15 s intervals and then at 
longer intervals until the strain became constant. 

Strain readings obtained on the two sides of the 
specimen were averaged and strain was plotted against 
time to obtain the loaded creep and unloaded recov- 
ery curves. Readings on the strain gages were nor- 
mally made to the nearest 2 X 10~ 5 in. Since a gage 
Jength of 0.25 in was used this is equivalent to a strain of 
8X 10" 5 . Thus, in the results given below, differences 
in strain of 1 X 1()" 4 are approximately equal to the 
minimum reading difference. 

Loads placed on the specimen (with a nominal 0.01 
in 2 cross-sectional area) varied from 5 to 40 lb giving 
stresses from approximately 500 to 4,000 psi. Most 
specimens were used for several runs, first at high 
and then at lower stresses. The first loaded creep run 
was considered a strain hardening treatment and data 
obtained on these runs were not used in the calculation 
of results other than for viscous strain rate. All runs 
were made at 23 ± 1 °C. 



3. Materials 

The specimens were prepared using a commer- 
cial alloy for dental amalgam certified to comply with 
American Dental Association Specification No. 1. 
This alloy (composition approximately Ag 70%, Sn 
26%, Cu 3.5%, Zn 0.5%) was mixed with mercury 
and condensed into a mold as described by Rodri- 
guez and Dickson [1] to produce a specimen with 
dimensions as shown in figure 1. Specimens were 
aged for at least one week to obtain essentially full 
mechanical strength [8]. 

4. Procedure 

The dumbbell-shaped specimen was placed in 
the grips and Tuckerman optical strain gages were 
mounted on opposite sides of the specimen as shown 
in figure 2. To obtain the creep curve, readings were 
taken on the strain gages, a weight was suspended 
from the lower grip and a second strain gage reading 




-0.08 n Rad. 



FIGURE 1. Dimensions of the dumbbell-shaped tensile specimen of 
amalgam. 



5. Results and Discussion 

The creep curves (both loaded and recovery curves) 
>f strain hardened dental amalgam as shown for a 




Figure 2. Tensile specimen in position for load application with 
optical strain gages mounted on opposite sides of the specimen. 
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number of different stress levels in figure 3 indicate 
that at room temperature amalgam exhibits three 
different types of viscoelastic phenomena: (1) in- 
stantaneous elastic strain, (2) retarded elastic strain 
(transient creep), and (3) viscous strain (steady-state 
creep). 

The viscous strain rate was determined from the 
loaded portion of the creep curve by taking the slope 
of the straight line portion of the curve, and was also 
determined from the recovery portion of the creep 
curve by dividing the value of the recovery strain 
(the permanent strain in the specimen) by the total 
time the load was on the specimen. The viscous strain 
rates for any given creep curve as calculated from the 
loaded and recovery portions of the curve were found 
to agree fairly well as shown in table 1. The log of the 
viscous strain rate was found to be a linear function 
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FIGURE 3. Creep and recovery of amalgam loaded in tension. 
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Straight lines are least-squares fits of equation log 6 r = log K + m log a to the data, with 
errors assumed to be in log e r only. Loaded values are from straight portion of creep curve; 
unloaded values are from strain remaining in specimen after unloading and recovery. 



of the log of the applied stress, as shown in figure 4; 
thus the viscous strain rate could be related to the 
applied stress by the following equation: 

e v = Ka»\ 12) 



€ v is the viscous strain rate, 
cr is the applied stress, and 
K and m are constants of the material. 

The value of m for amalgam is the value of the 
slope of the curve in figure 4, while the value of K 
is the antilog of the viscous strain rate value at a 
value of applied stress cr of 1 psi. For the dental 
amalgam used in this investigation values for K of 
2.85 XIO" 19 and 4.98 XlO" 19 were obtained from 
loaded and unloaded data respectively, and values of 
3.99 and 3.92 were obtained for m. 

The strain developed in amalgam due to the other 
two phenomena (1) instantaneous elastic strain and 
(2) retarded strain can be determined from the strain 
recovery since these two types of strain are recoverable 
while the viscous strain is not. Thus, at any given load 
the strain values taken from the creep curve after 
the sample has been unloaded (that is in the recovery 
portion of the creep curve) are subtracted from the 
strain value on the creep curve at the instant just 
before unloading of the specimen. This difference is 
plotted against recovery time U = Ti — T u where T u is 
the time at which the specimen was unloaded and 

Tarle 1. Viscous strain rates 



Specimen — run 


Stress 


Strain per minute 








lb/in 2 


Loaded Xl0- 5a 


Unloaded Xl0- 5h 


41-40-1 


3807 


5.65 


e 6.35 


41-40-2 


3807 


5.98 


5.80 


39-40-1 


3731 


4.92 


c 5.74 


39-40-2 


3731 


5.23 


5.08 


40-40-1 


3731 


7.45 


c 8.04 


43-35-1 


3458 


3.12 


c 3.71 


43-35-2 


3458 


3.23 


3 28 


43-35-3 


3458 


3.42 


3.44 


44-30-1 


2959 


1.92 


h 1.98 


44-30-2 


2959 


1.68 


1.68 


44-30-3 


2959 


1.74 


1.72 


41-30-3 


2855 


2.84 


2.96 


44-25-4 


2465 


0.91 


0.89 


35-25-1 


2375 


1.05 


c 1.33 


35-25-2 


2375 


0.98 


1.06 


35-25-3 


2375 


1.03 


1.0,8 


38-25-1 


2366 


1.09 


1 1 .40 


38-25-2 


2366 


0.783 


0.821 


38-25-3 


2366 


.785 


.806 


35-20-4 


1900 
1900 


.384 
.419 




35-20-5 


0.435 


38-20-4 


1892 


.349 


.358 


35-15-6 


1425 


.111 


.113 


38-15-5 


1419 


.076 


.077 


35-10-7 


950 


.031 


.034 


38-10-6 


946 


.023 


.025 


38-10-7 


946 


.017 


.018 


35-05-8 


475 


.008 


.007 



a From slope of straight portion of loaded creep curve. 
b From strain remaining in specimen after unloading and recovery. 

c These values include effects of strain hardening and were not used in calculating the 
relation between stress and viscous strain rate. 
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Tj is the time of the strain value on the recovery portion 
of the curve. These difference values, e\ are seen 
plotted against the recovery time for each load or 
stress in figure 5. These plots are a measure of the 
combination of the elastic and retarded elastic strain 
behavior of dental amalgam as a function of time for 
various stress levels. A measure of the combination 
of elastic and retarded strain may also be obtained 
from the loaded portion of the creep curve by taking 
values off the loaded creep curve and subtracting 
the viscous strain accumulated in the specimen at 
that time. The accumulated viscous strain at any time 
may be calculated by multiplying the viscous strain 
rate by the time corresponding to that value on the 
creep curve. Thus the difference between the creep 
curve value on the loaded portion and the viscous 
strain value at a corresponding time is a measure 
of the combination of the instantaneous and retarded 
elastic strain. However, a small error in the viscous 
strain rate causes a large error in the difference value. 
Therefore, the plot of the combination of elastic and 
retarded elastic strain versus time as obtained from 
the loaded creep curve is subject to large possible 
error. 

Using the approximation of Nakada's general 
formulation described by Leaderman, McCrackin, 
and Nakada [7], the ordinates of the inverted recovery 
curves as described above should be greater in the 
early portion of the time scale than the ordinates of 
the loaded portion of the creep curves minus the vis- 
cous strain. The difference between the curves should 
rise to a maximum within the first few minutes and 
gradually decrease to zero as time increases. However, 
as shown in figure 6, the loaded curves appear to lie 
above the recovery curves. Fitting such data to 



Nakada's formulation would be very difficult. More- 
over, the differences are not large and the method of 
obtaining values for the loaded curves is subject to 
considerable error because the loaded curves had to 
be corrected by subtracting the viscous component. 
Thus it is believed that the differences between the 
curves are not significant and so the curves were 
treated as though they are superimposable. Therefore, 
the data reported for the combination of elastic and 
retarded elastic strain were obtained from the recovery 
portion of the creep curves. 

The combination elastic strain (instantaneous and 
retarded elastic strain) becomes asymptotic with 
time in accordance with theory as seen in figure 5. 
The combination elastic strain values were plotted 
as a function of the various stress levels for cor- 
responding time, as shown in figure 7. The combina- 
tion elastic strain is seen to be a nonlinear function 
of the applied stress. When the combination strain 
values were divided by their corresponding stresses 
and then plotted against the corresponding stress 
for a fixed time, a linear plot was obtained for each 
fixed time as illustrated in figure 8; this result indicated 
that the combination elastic behavior of amalgam as a 
function of applied stress under the test conditions 
could be represented by an equation of the form: 



€' = A(t)(T + BHt)(TK 



(13) 



e' is the combination of elastic and retarded elastic 

strain, and 
a is the applied stress. 
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FIGURE 5. Strain recovery after release of tensile stress. 

Each curve is an average of 2 to 7 runs. 
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A(t) and B 2 (t) are functions- of time but not of 
stress. The value of A(t) for any time value is the 
intercept at cr=0 of the plot for that time value as 
shown in figure 8 while B' 1 {t) is the slope of the 
straight line for that time value. It is also noted in 
figure 8 that as a function of stress the combination 
strain divided by the stress is a straight line for all 
values of /. This indicated that over all ranges of t, 
the combination elastic strain obeys the same func- 
tional relation to the stress. Since the recovery curves 
were found to be superimposable upon the loading 
curves minus the viscous portion, it appeared that 
the nonlinear material did not obey the approximation 
of nonlinear generalization previously noted [7]. It 
was then questioned whether the retarded behavior 
of amalgam was truly nonlinear or whether the non- 
linear behavior could be attributed to the geometry 
of the specimen being observed. However, longer 
specimens were tested and found to give the same 
result, also photoelastic specimens were made and 
the stress distribution over the area observed was 
found to be uniform. It is therefore concluded that 
the observed nonlinear behavior of dental amalgam 
is not a geometric artifact, but an intrinsic phenomenon 
in the material. The difference between the recovery 
curves and the loading curves minus the viscous 
portion is given by [7): 

2o- ( r[0(A7\7 1 )-O)(7\7 , )|. 



where AT" is the time the specimen was under the 
load o- ()l . 
zero, so: 



load o" ,. For amalgam, this difference is found to be 



<&(A1\T) =<P(T.T) 



(14) 



for creep recovery from equilibrium, i.e., for A7" large. 
McCrackin [9] has suggested that this condition be 
satisfied by assuming that the value of <£> is only 
dependent on the value of the least of its arguments. 
That is, 



Table 2. Values oj 'A(t) and B 2 (t) in equation 

e' = M\)<r \ \\H\)(t- 



Time— min 


.-/(/) x 10 ? 


B*(t) x 10" 


0.25 


1.130 


1.906 


.50 


1.120 


2.468 


l.OO 


1.132 


3.065 


2.00 


1 . 1 83 


3.548 


3.00 


1.206 


3.876 


5.00 


1.248 


4.285 


10.00 


1.356 


4.737 


20.00 


1.518 


5.100 


30.00 


1.5% 


3.364 


40.00 


1.626 


5.636 


50.00 


L.648 


5.839 


60.00 


1.668 


5.970 


70.00 


1.693 


6.031 


80.00 


1.727 


6.050 


90.00 


1.761 


6.064 


120.00 


1.812 


6.223 


150.00 


1.877 


6.224 


2(X).0O 


1.952 


6.210 


300.00 


2.074 


6.082 


600.00 


2.103 


6.290 


looo. 00 


2.079 


6.497 



e' is combination of instantaneous and retarded 
elastic strain. 

<r is stress. 

Ait) and B*(t) are constants for any time and 
are functions of time but not of stress. 



The values for A(t) and B 2 (t) were determined 
by fitting curves to the data by the method of least 
squares and were tabulated as a function of time 
as shown in table 2. :i The 4{t) values were plotted 
as a function of the corresponding t values as shown 
in figure 9. The A(t) values are seen to approach an 
asymptote, A. ds as t— ► oo. Thus, the curve as shown in 
figure 9 could be represented by the following equation 
from linear viscoelastic theory since A(t) is the linear 
term in stress: 



AH 



~ J . + \?' 



(r)(l 



-Hr 



)dr, 



(18) 



<D(7',.r 2 ) = «-(r,)ifr, «r 2 . 



(15) 



With this condition, the quadratic term of eq (8) 
is shown in the Appendix to reduce to 
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BHT-0)dl<THd)] 



(16) 



so, eq (8) reduces to: 
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A(T-6)(1ct(6) 



\ 



<r(7') 



BHT-0)<1\(t-(»)\. 



(17) 



2 The notation B*(t) rather than Hit) is used for consistency 
notation of Leaderman. MeCrackin, and Nakada |7|. 
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Figure 9. Variation with time of the linear creep compliance 
term, A(t), in the equation e' = A(t)cr+ B 2 (t)cr 2 . 



•'Strain values for different stresses were obtained at each of the specific times listed 
in table 2 (and for many other times not listed in the table) by interpolation between the 
recorded values of strain along each of the 25 to 30 creep and recovery curves. Then the 
strain values for a specific time were fitted to the equation c'la = A{t) + B 2 (t)a to obtain 
the values of A(t) and B 2 (t) for that specific time. 
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XIO" 7 



A{t) =J + \J'-(T)dr- \ X J'{j)e-^dr, (19) 

J(k J 

A(t)=A as - f X ./'(r)e-"^T, (20) 



-A{t) 



J'(T)e-*dr, 



(21) 



where 

A(t) 

Jo 
J'(r) 



is the creep compliance term which is linear 

in stress, 

is instantaneous elastic compliance, 

is the retarded elastic compliance spectrum 

as a function of the retardation time r, and 

is the asymptote value of A(t) as t becomes 

very large. 



Thus, eq (20) does describe the behavior of A(t) as 
a function of time as seen in figure 9. The linear creep 
compliance term A(t) is plotted as a function of log t 
to obtain a signoidal curve as shown in figure 10. The 
first plateau of the signoidal plot at very small values 
of t should correspond to ,/ () (the instantaneous elastic 
compliance) [10]. The value of J calculated from the 
curve is 1.13 X 10~ 7 in 2 /lb, equivalent to a modulus of 
elasticity of 8.9 X 10 6 lb/in 2 . This value is in good 
agreement with a value of 9.1 X 10 H reported for the 
modulus of dental amalgam determined by an ultra- 
sonic method [11]. 

The nonlinear creep compliance term B 2 (t) was 
plotted against log t to obtain a sigmoidal curve as 
shown in figure 11. Thus, in the nonlinear theory of 
Nakada [6] the experimental B(t) for amalgam might 
be represented by the following equation: 



Bit) 



= P/3(y) 

Jo 



(I" 



,-tly) 



dy (22) 



which would indeed describe the behavior of the curves 
seen in figure 11. The combination elastic behavior 
of dental amalgam in creep under the test conditions 
used could then be described by means of the follow- 
ing equation from viscoelastic theory [6, 10]: 

e'=A(t)(T + B 2 (t)o-\ 



-J or + <r 



J f (T)(l-e-^)dr 



+ 



f X /3(y) 

Jo 



(i- 



,-tly 



)dy 



(24) 



and since the loaded portion of the creep curve for 
amalgam is also composed of viscous strain e v as well 
as combination elastic strain e' then the strain on 
the loaded portion is composed of the sum of the two 
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FIGURE 10. Variation with log of time of the linear creep com- 
pliance term, A(t), in the equation e' = A(t)cr+ B 2 (t)cr 2 . 
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FIGURE 11. Variation with log of time of nonlinear creep compliance 
term, B 2 (t), in the equation e' = A(t)(j- + B 2 (t)cr 2 . 
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folic 



= €' + €„. 



(25) 



+ *±k(*=$*) w 



(27) 



Thus 

e = ./ cr-ho- J J'{r){\-e-^)dr 



+ 



i; 



/8(y)(l-e-^)rfy o- 2 + /Co-'"* 



To- 2 - 



(26) 



Therefore, applying both linear and nonlinear visco- 
elastic theory to the experimental behavior of amalgam 
under the test conditions used, a general equation 
describing the creep behavior of amalgam as given 
by eq (26) above is obtained. 

Applying the nonlinear generalization of the Boltz- 
mann superposition principle as developed by Nakada 
[6], to the above creep equation for amalgam, the 
stress-strain curves for various stress rate conditions 
were calculated for amalgam from the creep data and 
compared to the experimental stress-strain curves 
obtained under those conditions. The stress-strain 
curve corresponding to the experimental stress rate 
was calculated from the following equation: 



E,(nAt) = ^ (o-,-o-/-,) 



4/1-1+^,, + i 



<) 



+ 2 (°"?-°f-i) 



Bn-\ +#,7+i-i 



which is a finite difference approximation to eqs (17) 
and (12). For the experimental data, the testing 
machine was run at a constant head speed which 
produced a varying stress rate over a total time of 
3V2 min. From the stress versus time curve, observed 
stress values were obtained at 15-s intervals. Using 
the values of A(t) and B(t) calculated for 15-s intervals 
from the creep data and also the viscous constants 
from these data, a curve was calculated for comparison 
with the experimental curve. Figure 12 shows the 
contribution of instantaneous elastic, retarded elastic 
and viscous components to the total calculated strain. 
It is seen in figure 13 that good agreement is obtained 
between the calculated and experimental stress-strain 
curves. Thus, it is concluded that in the case of dental 
amalgam the results of creep tests can be related to 
those of stress-strain tests by use of viscoelastic 
theory [6, 10]. 



6. Conclusion 

Dental amalgam was found to exhibit three types of 
viscoelastic phenomena: (1) instantaneous elastic 
strain, (2) retarded elastic strain (transient creep), 
and (3) viscous strain (steady state creep). 

The combination of elastic plus retarded strain 
in tension was found to be a nonlinear function of 
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FIGURE 12. Calculated Instantaneous elastic strain E l5 retarded elastic strain E R , 
viscous strain E v , and total strain E T occurring during a stress-strain test of dental 
amalgam. 

Head speed of testing machine was constant, and total time was 3V2 min. 
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stress and could be represented by an equation of 
the following form: 

e' = Act + B 2 cr 2 where A and B are functions of time, 

but not of stress. 

The viscous strain rate is also a nonlinear function 
of stress and can be represented by an equation of the 
form: 

€v = Kcr m where K and m are constants of the material. 

An equation describing the creep behavior of 
amalgam was obtained by the application of nonlinear 
viscoelastic theory to the experimental behavior of 
amalgam under the test conditions used. By applying 
the nonlinear generalization of the Boltzmann super- 
position principle to this equation the results of 
creep tests were related to results of stress-strain tests. 



7. Appendix 

The derivation of eq (16) is as follows [9]: 

r cr(T) 



The first integral is 



fcr(T) r Caid) "1 

B 2 (T-0) dxr(<l>)\<kr(0) 



CaiT) 

-J. ' 



B 2 (T-6)o-(6)d(r(0) 



1 fcr(T) 



B\T - 6)d[& 2 {d)\ 



Applying Dirichlet's formula [12] to the second 
integral, 



r \r rut ,n ^)^) ,f 



d<j} dO 



dO 



gives 

r r[V(r-4>) 



da(<f)) rfcr(fl) 
d<f) dd 



d0\d4> 



CaiT) Cat 
Jo J() 



<5>(T-d,T-<t))d(r(0)d(T((t)) 



CaiT) Ca{6) 
Jo Jo 



B 2 (T-d)da((f))d(T(6) 



+ 



Ca{T) CaiT) 

1 

Jo J a{0) 



B 2 (T-(f))d(T(<t>)d(T(0). 



do-(cf)) f*d<r(6) 



fh^^j,,^ 



P 

Jo 



d(f> 



f T B 2 (T-<t>) 

Jo 



da(<f>) 
d<f) 



cr((/>) dcf) 



-*£»*-*■>% 



2 (</>)] 



dcf) 



d(f). 
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FIGURE 13. Comparison of calculated and experimental stress-strain curves. 
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However, this is the same as the first integral as 
6 and are dummy variables of integration so the sum 
of the integrals is 



Jo 



(T-e) 



d[<r 2 (6)] 



10- 



f<r(T) 

■J. "' J 



do 
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